Dynamic Response of Adhesion Complexes: Beyond the Single-Path Picture 
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We analyze the response of molecular adhesion complexes to increasing pulling forces (dynamic 
force spectroscopy) when dissociation can occur along either one of two alternative trajectories 
in the underlying multidimensional energy landscape. A great diversity of behaviors (e.g. non- 
monotonicity) is found for the unbinding force and time as a function of the rate at which the 
pulling force is increased. We highlight an intrinsic difficulty in unambiguously determining the 
features of the energy landscape from single-molecule pulling experiments. We also suggest a class 
of "harpoon" stickers that bind easily but resist strong pulling efficiently. 

PACS numbers: 82.37.-j, 87.15.-v, 82.20.Kh, 33.15.Fm 



The last decades have witnessed a remarkable devel- 
opment of physical investigation methods to probe sin- 
gle molecules or complexes by various micromanipula- 
tion means. New techniques have been put forward 
to probe the unfolding of proteins and to quantify the 
strength of adhesion structures . An important step 
in this direction is the proposal of the group of Evans 
to use soft structures to pull on adhesion complexes or 
molecules at various loading rates (dynamic force spec- 
troscopy) j|] . Moving the other end of the soft structure 
at constant velocity induces on the complex a pulling 
force that increases linearly in time f = rt. Measur- 
ing the typical rupture time ttyp yields a typical rupture 
force /typ = rttyp that depends on the pulling rate r. 
This provides information as to the energy landscape of 
the bound complex. Indeed, in many situations one ob- 
serves a linear increase of /typ with log(r), which can be 
understood within a simple adiabatic Kramers picture 
for the escape from a well (bound/attached state) over 
a barrier of height E located at a projected distance x 
from the well along the pulling direction. The progressive 
increase of the force results in a corresponding increase 
of the escape rate, so that, in agreement with some ex- 
periments the typical rupture force increases loga- 
rithmically with r: /typ ^ /cbT/x ln[rx/(/cBTcc;)], where 
LO is the escape rate in the absence of force. The rupture 
time on the other hand decreases with r. The occur- 
rence in some cases of two successive straight lines in a 
[/typ, log(r)] plot has been argued to be the consequence 
of having two successive barriers along the ID escape 
path, the intermediate one showing up in the response 
at fast pulling rates [|j (Figs, [^a and|2|). Other theories 
have tried to back up more complete information as to 
the overall effective ID potential landscape by an analy- 
sis of the probability distribution for rupture time and of 
the statistics of trajectories before rupture [|^,|[. Assem- 
blies in series and in parallel of such ID bonds have also 
been considered ]9|-pT[. 

In this Letter we point out limitations arising from 



the a priori assumption of a single-path topology of the 
energy landscape for the interpretation of such exper- 
iments. From the analysis of simple examples with a 
two-path topology, we draw three conclusions: (i) first, 
the dependence of the rupture force and rupture time on 
the pulling rate can take various forms, including non- 
monotonic behavior (see e.g. Figs. || to ||). (ii) Second, 
the main features of the energy landscape can not be un- 
ambiguously deduced from a [/typ, log(r)] plot, as very 
different landscapes can yield similar curves (Fig. o). (iii) 
Third, we propose simple "harpoon" designs (Fig.]l| c and 
d) for functionally efficient stickers that can bind easily 
but resist strongly in a range of pulling forces (Fig. 
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FIG. 1. Sketch of the topology of the main valley of the 
energy landscape for a few examples. denotes the fundamen- 
tal bound state, A and B are local minima, and a, a , b and 
b' are passes to overcome. To the right (increasing values of 
x) of the last passes is the continuum that describes unbound 
states, (a) classical single-path scheme. (b,c,d) unbinding can 
occur through two alternative routes a and p. 
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FIG. 2. Classical picture for a single-path energy land- 
scape (Fig. 0a): the probability density P{f) for unbinding 
at force / is plotted in grey-scale as a function of the pulling 
rate r. The typical force /typ (locus of the maximum of P) 
is highlighted with a dashed-line. Plotted curves correspond 
to 12, x«/ 0.5, Ea = 9, xa = 1, and Ea 20, 

Xa = 2. At very low pulling rates unbinding is not affected by 
the pulling and proceeds over barrier a with a "spontaneous" 
rate a;oexp(— ^a)- For larger pulling rates the typical un- 
binding force /typ increases linearly with log(r), with a slope 
proportional to 1/xa- Increasing further the pulling rate can 
lead to a steeper slope oc 1/xa' corresponding to escape over 
the inner barrier a . These asymptotes are depicted with solid 
lines. The dashed arrows along the drawings indicate which 
pairs of energy well and barrier are probed in these asymptotic 
limits. Inset: mean rupture time against pulling rate. 

Obviously for real binding/ adhesion complexes, there 
are numerous (conformational) degrees of freedom, and 
the configurational space is clearly multidimensional. 
This allows for complex energy landscapes and vari- 
ous topologies for the structure of their valleys and 
passes [|l2[. Only the probing (pulling) is unidirectional. 
We note in passing that even for more macroscopic sticky 
systems, usual adhesion tests for soft adhesives often 
show up hysteresis loops associated with the existence of 
more than one degree of freedom jl^] . We do not attempt 
here an exhaustive exploration of effects allowed by the 
multidimensionality of the phase space, but rather focus 
on a few simple two-path topologies (Fig. [l]), to argue for 
the three points mentioned above. 

The three examples we consider, sketched in Figure 
|l| b, c, and d, correspond to simple hairpin schemes 
whereby detachment can proceed through two alterna- 
tive routes a and f3. These simple quasi ID schematic 
situations can be conveniently dealt with using an adi- 
abatic Kramers theory, which has been shown to be an 
efficient way of obtaining semi-quantitatively correct an- 
swers 

A common set of notations can be ascribed for all cases 
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FIG. 3. Switch geometry (Fig. |lp): plot of the same quan- 
tities as in Figure H, for Ea — 20, Xa — 0.5, and E^ — 30, 
Xb = 2. At low pulling rates unbinding is controlled by es- 
capes over a whereas for large values of r it occurs mostly 
over b: the slope of the unbinding force (average or typical) 
decreases from 1/xa to 1/xh. 

(Fig. 0). From the fundamental bound state "0", the 
route a for escape (detachment) is over barriers a, of 
height Ea located at a projected distance Xa from "0". 
Alternatively, escape can occur through branch over 
barrier 6, of height £^5 and projected distance X5. All 
energies and projected distances are measured relative 
to the state "0" (i.e. £^0 = and xq = 0). Inter- 
mediate barriers a', b' and local minima A and B may 
exist, with energies Ea' ^Ei)' ^Ea-,Eb (all positive), and 
projected distances x^', x^, ^b- In line with typical 
values from experiments, we choose to write energies in 
units of /cbT 4 pNnm and distances in nm. 

Practically, we describe the time evolution of the prob- 
abilities of being in the potential minima (bound states) 
using "chemical" transition rates over the barriers as 
given by the Kramers formula. We furthermore assume 
the attempt frequencies to be constant and all equal to uoq 
which provides the only intrinsic time-scale in the prob- 
lem, so that the transition rate from minimum / over the 
neighboring barrier i is uoq exp[— (£—£/) + /(t)(xi —x/)]. 
For the plots of Figures ^ to || we take arbitrarily ujq = 
10^ s~^. Jump over the rightmost barrier (a or b) of 
either path corresponds to rupture leading to escape to 
X ^ 00. 

We focus on the case where either £5/ or £5 is larger 
than Ea so that a is the "natural" route by which attach- 
ment and detachment proceeds in the absence of pulling. 
We also limit ourselves to simple scenarios in which the 
force is linearly increased in time / = rt. 

For further reference we recall the classical single-path 
scenario (Fig. |^a) in Figure || for a typical set of param- 
eters, and then we turn to a brief analysis of the three 
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FIG. 4. "Harpoon" geometry (Fig. |l|c): plot of the same 
quantities as in Figure bl for Ea = 20, Xa = —2, and Eb = 40, 
= 2. Pulling here impedes unbinding through the "spon- 
taneous" route a, so that as soon as the rate is strong enough 
for pulling to affect unbinding, the escape is controlled by the 
larger barrier 6, resulting in an upward jump of the typical 
unbinding force and time. Inset: the average unbinding time 
is here non-monotonic. 

geometries we have introduced (Fig. |l| b, c, and d). 

First case: switch - Topology as in Figure |l|b. Escape 
occurs through either barrier a or barrier b both located 
downwards in the pulling direction {xa^Xb > 0). The 
escape proceeds through path a at weak pulling rates as 
Ea < £^6, but if > Xa it can switch to path P for pulling 
forces / large enough such that Ea — fxa > £5 — /x^. The 
result (see Fig. ||) is then a succession of two straight lines 
of decreasing slopes in the [/typ, log(r)] plot, the first one 
(slope oc 1/xa) characteristic of the spontaneous route 
a while the second (slope oc I/X5) provides information 
on the alternative route p. In the trivial case Xa > 
route P is never explored so that the classical single-path 
picture applies. 

To clarify the calculation leading to the plot in Figure 
H, we describe the evolution of the probability of attach- 
ment p{t) at time t of the system initially attached at 
time t = [p{0) = 1] by 

dtp{t) = -ujoie-^-^^^'^''- + e-^^+^(')^^)p(t) (1) 

Solving (1) numerically with f = rt yields p(t) and there- 
fore the probability density P{f) = —\dtp{^) for the 
unbinding force. The typical values of / are highlighted 
in the plots, with the whole distribution P{f) suggested 
through a grey-scale. Similar procedures will be used 
in the following examples, with thermal equilibrium be- 
tween the bound states assumed as initial conditions. 

Second case: harpoon - Topology similar to the pre- 
vious one but with < (Fig. |c). The main feature 
here is that as the pulling force increases, the probabil- 
ity to escape over a decreases. Therefore the system gets 
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FIG. 5. "Selective harpoon" from the combo topology 
(Fig. ^): same quantities as in Figure ^, for Ea — 20, Xa — 2, 
Eb' 10, xy 0.5, Eb = 5, = 1.5, and Eb 27, Xb 2.5. 
At low pulling rates the spontaneous path a is used. Upon 
increase of r, larger forces are employed and the minimum 
B becomes favorable as compared to 0. As Ey is not too 
large, equilibration of population then empties in so that 
escape eventually occurs from B over 6, resulting in a higher 
straight line of slope oc l/{xb — xb)- At even higher pulling 
rates, because Xa > xy ^ the escape over a becomes faster 
than this equilibration, and therefore, path a is used again. 
Barrier a controls the behavior at low and high rates, but 
in an intermediate window, a stronger bonding is provided 
by barrier b. The typical (dashed line) or average unbinding 
force is non-monotonic. 

"stuck" in route p. If the barrier £^5 is infinite (left side of 
Fig. |l|c), there is a finite probability Poo = exp(— ) 
that unbinding never occurs. For a finite but high barrier 
£5, pulling eventually results in unbinding but at high 
rupture forces (see Fig. ^). The topology thus allows 
here to form "easily" (i.e. over barrier a) a "harpoon" 
sticker that can resist strong pulling. Correspondingly 
the mean unbinding time increases first with pulling rate 
(a phenomenology connected to the negative resistance 
analyzed in Ref . ) , before decreasing for larger values 
when activated escape over b dominates. Note that the 
probability distribution P{f)^ now consists of two sep- 
arate ensembles, which coexist over a narrow region of 
pulling rates. This is in contrast with Figure ^ where 
there is a continuous evolution of a single cloud. 

Third case: combo - The alternative route consists of 
two barriers and a local minimum B (Fig. |l|d), and we 
focus on the case where £5/ is smaller than the two others. 

Thanks to the increased complexity and number of pa- 
rameters in this case many scenarios can occur, covering 
features already unveiled in Figures § to || (e.g. switch 
and harpoon). More intricate pictures can also show up, 
as depicted in Figure |[ An explanation of this example 
is given in the caption, illuminating how for low or high 
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data of Figure |2|. (b) classical single-path of Figure 
Ea' — 11, Xa' — 1, Ea = 8, xa = 1-5, and Ea — 20, Xa 
(c) combo two-path geometry (Fig. |l|d), with Ea - 
Xa — 2, Ey — 18, xy — 2.5, Eb — 15, xb — 3, and Eh - 
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FIG. 6. Similar curves obtained from significantly differ- 
ent energy landscapes, (a) classical single-path of Figure |l|a, 

with 
= 2. 

20, 
= 27, 

Xb = 3.5. In all cases the straight part for weak r corresponds 
to escape over a from 0. The second steeper slope corresponds 
to escape over a from in case (a), over a from A in case 
(b), over b from B in case (c). Neither the topology, nor the 
location of the probed segment of the energy landscape can 
be asserted from such data sets. 

pulling rates barrier a controls the behavior, whereas for 
intermediate values, the secondary and stronger barrier 
h limits unbinding. Two features are striking. First, the 
unbinding force (typical or average) is no more mono- 
tonic. Second, branch 13 results in a strengthening of the 
binding complex for a given window of pulling rates r 
(selective harpoon). 

Discussion -WiXh the three simple examples above, we 
have clearly enlarged the numbers of behaviors that one 
may obtain from a classical dynamic force spectroscopy 
method (see Figs. || to ||). Conversely, we also want to 
stress the second point (ii) mentioned in the introduc- 
tion: simple patterns (e.g. the succession of two lines of 
increasing slopes) can be the outcome of many diverse 
landscapes. For example. Figure ^ displays force-rate 
curves similar to that of Figure ||, but that correspond 
to sensibly different landscapes. Not only are the typical 
and average unbinding forces very similar in the three 
cases, but so are the probability distributions for most 



values of r. Only close to the cross-over between the two 
straight lines can slight differences be detected. To dis- 
tinguish more selectively possible landscapes, it may be 
necessary to use other temporal sequences than the sim- 
ple / = rt, e.g. to reveal equilibration processes between 
local minima. 

Eventually we would like to emphasize that the har- 
poon geometries proposed here constitute a very obvious 
paradigm for efficient stickers. Attachment of the sticker 
can proceed through route a with a possibly not too high 
barrier Ea- The harpoon configuration then allows to 
benefit from the much stronger h barrier for a given win- 
dow of pulling forces, making the sticker more efficient 
in these conditions. This "hook" design is obviously a 
favorable strategy for adhesion complexes, the function 
of which is to maintain adhesion under the action of well- 
defined tearing stresses. It would be surprising if advan- 
tage was not taken of this by some biological systems. 
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